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^3 ■ 1 Introduction 

Professor Erber has made important contributions to several areas of both pure and applied physics, making 
it easy to identify topics about which one can write to celebrate his work; approaching such topics at his 
level of rigor and insight is far more difficult! His contributions to fundamental electromagnetic theory and 
q-i quantum mechanics in particular include papers on electromagnetic energy density in dispersive media [T], 

synchrotron-Cerenkov radiation [2], radiation reaction [3], and quantum jumps 0J, to cite only a few of those 
with which we are familiar. The first two papers cited, among others, deal with electromagnetic processes 
■ in dielectric media, and a small part of that subject will be addressed here. Specifically, this article is 

Q\ | concerned with the quantized electromagnetic field in a dispersive and dissipative dielectric medium, and 

with the energy density in such a medium. Erber's work has also dealt with nonlinear processes in such 
media as well as in vacuum, but we will restrict ourselves here to linear and idealized, homogeneous media. 
This is an important subject about which much has been written, although in most of the literature it 
(N ■ has been assumed that the medium can be assumed to be non-dissipative at field frequencies of interest. The 

theory ignoring dissipation is not without value; it can be used to describe, for example, the spontaneous 
emission by an atom in a host dielectric that does not absorb radiation at the atom's transition frequency. 
But a dispersive medium cannot be non-absorbing at all frequencies. What happens, for instance, if the 
medium in our example strongly absorbs radiation at the transition frequency? And what about situations in 
which we cannot exclude any frequencies a priori and therefore cannot ignore absorption? In the calculation 
of the van der Waals force between two neutral dielectric bodies, for example, all field frequencies can in 
principle contribute to the force. For two perfectly conducting parallel plates at zero temperature, similarly, 
Casimir [5] discovered, as a consequence of the zero-point electromagnetic energy of every mode, that there 
is an attractive force per unit area between the plates. His original method involving changes in zero-point 
field energy was later extended to dielectrics by van Kampen et al. [6] and others [7]. As emphasized by 
Ginzburg [8J, however, these theories invoking changes in zero-point energy ignore absorption: "... oddly 
enough there is no mention that they consider directly only transparent media" [9|. An entirely different 
route, based on the calculation of the force via the stress tensor, was taken by Lifshitz [TU]; he accounts for 
absorption through the fluctuation-dissipation relation between the quantum fluctuations of the polarization 
density and the imaginary (absorptive) part of the permittivity. 

The intent here is to derive, in probably the simplest way possible, expressions for the quantum elec- 
tromagnetic field in a dispersive and dissipative dielectric medium, treating the medium as a continuum. 
The derivation might seem superfluous in the sense that correct expressions for the electric and magnetic 
fields in such a medium are already available [llj . However, the diagonalization procedure by which these 
expressions are obtained is not easily applied to general, inhomogeneous media, whereas the "Langevin ap- 
proach" presented here can be applied more or less straightforwardly when extended and formulated via 
Green functions [T2l. 
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We begin in the following section with the simple model of an oscillator A coupled to a reservoir R of 
other oscillators, the R oscillators causing a damping of A described by a Langevin equation. In Section [3] 
we review the Fano procedure for the diagonalization of the Hamiltonian of this system, and compare the 
diagonalization and Langevin-equation approaches insofar as they describe the time evolution of A. In the 
limit of zero temperature, or at any finite temperature, the two descriptions are shown to be equivalent. In 
Section[5]we generalize these considerations, following the Langevin approach, to the model of a homogeneous 
dielectric medium in which each atom is treated as a harmonic oscillator coupled to its own reservoir. The 
Langevin forces acting on the atoms give rise to a noise polarization determined by the reservoir operators, 
and the fluctuating electromagnetic field caused by this noise polarization can therefore be expressed in terms 
these operators. In similar fashion to the model of Sections [5] and O the quantized electric and magnetic 
fields obtained in this way have exactly the same form as obtained by Fano diagonalization. It is shown 
explicitly in Section [5] that the zero-point energy per mode of frequency ui is (l/2)hu> regardless of whether 
there is absorption at that frequency [13] . 



2 An Oscillator and a Reservoir 

Consider an oscillator A of frequency ujo coupled to a reservoir R of other oscillators, a well-studied model for 
dissipation in quantum theory. To make things as simple as possible we will assume that the A-R coupling 
involves only energy-conserving processes, and choose this coupling such that the Hamiltonian is 

/>oo />oo 

H = fiui^a) a + / du>hujw (ui)b(uj) + hy^j/n / dui[a)b(uj) + b\u>)a], (1) 
Jo Jo 

with [a, e^] = 1, [b(u>), b(u>')] = 0, b'(oj')] = S(u> — u'). (We use the circumflex to denote operators.) 

Because only energy-conserving processes are included, our model is consistent with the so-called "rotating- 
wave approximation" (RWA) . The coupling we have chosen results in a frictional damping rate 7 in the time 
evolution of A, as follows. 

The Heisenberg equations of motion for a and b(u>) are 



a = —iujoa — iy/^f /n / diob(uj), (2) 
Jo 

b(uj) — —iuib(uj) — iy/'y/na. (3) 

Using the formal solution 

b(w,t) = h{oj)e- iut - f dt'a{t')e iui< t'-V (4) 

Jo 

of equation ([3]) in equation and defining &o(w) = o(w,0), we obtain 

h{t) + iu a(t) + ( 7 /tt) / duj \ dt 'a{t')e luj{t = -iyfiyfc I dujb (u)e- iuit , (5) 



the operator on the right-hand side being a quantum Langevin force. We solve this equation for "steady 
state" ("ft 3> 1) by first writing 

poo 

a(t) = / dnA{n)b Q (n) e - int . (6) 

Jo 

Then, using the approximation 

* , i 

tfJbo-W -t) s n g, u _ j« _ ip /fl 

„ LO — \l 
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for times t such that fit ^> 1 for frequencies Q that make a significant contribution to the time evolution of 
a(t), we obtain 

A ^ = n-J^n) + i 7 - (8) 

The frequency shift 

m =i P r-^- ( 9) 

obviously diverges in our model, and so the upper limit of integration must be appropriately cut off; for our 
purposes there is no need to explicitly indicate any cutoff. Then our solution for a(t) is 

W V IT J n-LU + A(fi) + i7 v ' 

Note that, under the assumption that ojq 7, consistent with the RWA, we have 



[&(t), a T (t)] = -1 / ? ^— - = 1, (11) 

7T J Q (U- U>qY + Y 



as required for the validity of the RWA. 



3 Fano Diagonalization 

A method of diagonalizing a Hamiltonian for coupled oscillators, used many years ago by Fano [14] . has been 
employed in seminal work by Huttner and Barnett [11] to obtain expressions for the quantized fields in a 
dissipative dielectric medium. Here we briefly review the method in the case of the model Hamiltonian |1| , 
and compare it to the approach of the preceding section. 
We define an operator 



B(Q) = a(0)a + / dw^(0, u)b(u) (12) 
Jo 

that we require to satisfy 

[B(fi),&(ci')]=s(n-n'), [B(n),B(n')} = o, (13) 

and also require that the Hamiltonian (Q} takes the diagonal form 

P 00 

h= dmn&(n)B(n). (14) 



From [B(fl),H] = h£lB(n) and the definition (fT2")l we deduce equations relating the coefficients a(fi) and 
/J(n,w): 

(w - ft)a(fi) = J dLuP{Q,Lu), (15) 

(cj - fl)/3(Cl, u) + J^(fl) b(ui) = 0. (16) 



The last equation leads us to write 



and it follows from (11511 that 



>'(0^i ( .(())/(.<)).>V ^ - A/Z^Z^: ( 17 ) 



/(n) = J^[n-wb + A(n)], (is) 
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i.e., _ 

P(Sl,u) = ,&(fl)[fl - Ub + A(fl)]% - fl) - •/I^, (19) 

y 7 v k & — it 

where A(Q) is defined by ([9]). 

To determine a(Sl) we impose the requirement that the commutation relations fp~3|) be satisfied. From 
the commutation relations stated earlier for the a and b operators we obtain 

OO />OC 



[B(Sl),Bl(SY)] = a(Sl)a*(Sl')+ / du dw'/3(n, u)/3*(Q', u/)5(u - (J) (20) 

Jo Jo 

or, from (fl9]) and some straightforward algebra, 

[B{n),B\n')} = - a (si) a *(si') x j[si-w + A(si)] 2 <s(si-si') 

t 2M + VrA P r^- % -„<)}. (2i, 

7T Si'— il 7T Z Jq CJ — Si Jo U) — il J 

Using 

P—^ = k + iirS(uj — Si) (e^0+) (22) 

and partial fractions we obtain 

\B(n),BUn')} = -a(fi)a*(n') (7 2 + [fl — u>o + A(S1)] 2 ) <5($7 — Si'). (23) 
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Therefore we can satisfy (|T3|) by taking 



Then 



a ^ = « — Afo^ ■ ■ (24) 

\l — ujo + A(S2) — 17 

a(t) = [°° dQa*(Q)B(Q,t) = ./l /°°dO , , (25) 

Jo V I" Jo SI - w + A(S1) + 

since B(Sl,t) = B(Q, 0)e-* nt . 

The expressions (|10l) and (|25[) obtained respectively in the Langevin and Fano approaches look formally 
the same in the sense that [B(Q, 0), B+(Sl', 0)] = [6 (Sl), SJ(sT)] = <5(Sl - Si') and [5(S1, 0), B(Q\ 0)] = 
[6o(M), 6o(^')] = 0- They differ in that B(fl, 0) in (|25[) is a linear combination of A and R operators, whereas 
only R operators determine a(t) in (jTTJjl . Suppose, however, that at t = the A and R oscillators are all in 
their ground states. This state |\&) is the exact ground state of the coupled A-R system in the RWA: 

fr(n,t)B(n,t)\y) = 0. (26) 

In this case the properties of A derived from (fTU|) are trivially equivalent to those obtained from (f2"5|) . If 
the system is not initially in an eigenstate of B^B, it will nevertheless approach after a time 3> 7 -1 an 
equilibrium state for which the long-term solution (fT0|) for a(t) is applicable, i.e., transient effects associated 
with a(0) at some initial time t = ultimately play no role in the evolution of A. Equilibrium values of 
correlation functions involving products of the B operators are determined solely by the reservoir operators; 
in thermal equilibrium, for example, (B^(Cl, t)B(Q', t)) = (b^(fl,t)b{Cl',t)) = [exp(7iSl/fc s T) - l]- 1 ( 5(Sl - Si') 
and it follows from either (p~0|) or (|25[) that (aJ(t)a(t)) — [exp(hujo/kBT) — when we invoke the condition 
luq 3> 7 as in (jllj) . In other words, in any state of equilibrium the solutions (|10[) and (|25[) provide equivalent 
descriptions of A. This equivalence holds more generally beyond the RWA (see below). 
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4 Fields in a Dielectric Continuum 



Aside from the need to introduce oscillator strengths in order to obtain correct numerical results, we can 
model a dielectric medium in which atoms remain with high probability in their ground states as a collection 
of harmonic oscillators. We take each oscillator "atom" to have a mass m and a natural frequency uiq. We 
assume each of these material oscillators is coupled to a reservoir of other harmonic oscillators responsible 
for the damping of its oscillations and line broadening. For the Hamiltonian, including the electromagnetic 
field and its (electric-dipolc) coupling to the material oscillators, we write 

poo ■ 2" r°° 

+ duTwJ2 bt(w)-b i («) + - -i / dwA(w)^% • [bj(cj) - bj(cj)]. (27) 
Jo . L 2j J . 

The first two terms are the Hamiltonian operators for the electromagnetic field, the material oscillators 
(atoms), and their coupling via the (Coulomb-gauge) vector potential A(rj), rj being the position of the jth 
atom. The third and fourth terms are respectively the Hamiltonian operators for the reservoir oscillators 
and their interaction with the atoms. The reservoir oscillators satisfy the commutation relations 

[b ifl (u>), &}>')] = SijS^Siw - «'), [6i„(w), b jv (w')] = 0, (28) 

where we use Greek letters to denote Cartesian components of vectors. We choose the atom-reservoir coupling 
constant to be 

mh'juj \ 



A(u>) = I— J- J . (29) 

Then, as shown below, each atom's oscillations are damped at the rate 7. Note that no rotating- wave 
approximation is made in writing (|27[) . The operators A and E satisfy the usual canonical commutation 
relations for the electromagnetic field. 

From (J25J) and [ iipipju] — ifaSijSfiv we obtain the Heisenberg equations of motion 

L+u&iij = — E(r») + — / dtA(w)[b,-(w,t) - bt(w,t)], (30) 
m m J J 

bj(o;,t) = — icjbj(w,t) + — A(w)xj . (31) 

Using the formal solution of (|3Tj) in (|30| . we write 

e - 1 - if 00 /"* 

= — E(rj) + — F Lj (t) H diuA 2 (cu) dt'xj(t')[2i sin w(t' - t)], (32) 

to 777 mn. Jo Jo 

where the Langevin force operator F Lj (t) acting on the Jth atom is 

/>oo 

F Lj (t) = i dcuA{uj)[b 3 (uj, 0)e- iut -b] (w,0)e iart ]. (33) 
Jo 

The third term on the right-hand side of (132)) is 



(2wA 2 (u) / dt'x^t') sia a; (i' - t) = / dt'x.j(t') 



rah Jo 

du)u)sinw(t' -t) = 2j dt'xj(t')—S(t' - 1) = --yxj(t) . (34) 

Jn Ot 
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We ignore a divergent frequency shift which, as in the model considered in the preceding sections, can be 
made finite by introducing a form factor or a high-frequency cutoff; the (finite) shift can be assumed to be 
contained in the definition of ojq. Equation (|32p then has the form of a quantum Langevin equation: 



,2,-, 



1 



-r 7x4 + wjx,- = — E(r,-) + —F L j(t). 

m m 

In the absence of coupling to the electromagnetic field we have, for times t ^> 7 -1 , 

bUuj)e luJt 



Pj(t) — mxj(t) = / duojK (uj) 



bj{u)e- 



uJq — lo 2 — ijuj uj 2 — lo 2 4- iju 



(35) 



(36) 



(We now write bj(u) in place of bj(aj, 0).) Similarly, using (|28[). we obtain 



[Xiii (t),Vjv{t')\ 



duiui 2 cos uj(t' — t) 

(lUq - UJ 2 ) 2 + J 2 U! 2 



COSCt>i(i' — t) 



7 

2wi 



sinwilt' — <| 



-t|*'— t|/2 



(37) 



where u/i = [uJq — 7 2 /4] 1 / 2 . The canonical commutation relation [ij /i (t),p J v(i)] = ihSijS^u is therefore 
preserved in the coupling of each atom to its reservoir. 

Since we are working in the Heisenberg picture, expectation values are over the initial state of the 
coupled system of oscillators. If we assume that the reservoir is in an initial state of thermal equilibrium at 
temperature T, then 



(&]>)& 3 v(w , )> = (M^)& >')> - M^(«> - 



f 



SijS^S^u - to'). 



(38) 



The Heisenberg equations of motion for the electric and magnetic fields that follow from the Hamiltonian 
(|27| and the canonical commutation relations for the field operators have exactly the same form as their 
classical (Maxwell) counterparts: 



V x E 



ldB 

"cdt' 



47T - 

V x H = — J 

c 



1 9E 

c~dt' 



(39) 



For a charge-free medium, furthermore, V • B = V • D = 0, where 

D = E + 4ttP, 

0P(r,t) 



J(r,t) 



at 



P(r,<)= e ^x J (^ 3 (r-r J ), 



(40) 



with B = H in our model. 

It is convenient to work in the frequency domain and write 



/>oo 

E(r,i) = / duj[E(r,oj)e- luJt +E t (r,w)e lu 
Jo 

/>oo 

H(r,i) = / rfw[H(r,w)e-^ t + H t (r,w)e l 
j 



P(r,t) = 



du[P(r, u>)e 



■P + ( 



r, w)e 



(41) 
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The Fourier transform of the polarization density may be written as 

P(t,w) = e^x.H^r-r,), (42) 



*j(t) = / du[x j (u)e- iut +i&(u)e iut ], (43) 



and it follows from (1351) that 



P(r >W ) = e V m VE(r J ,^ 3 (r-r J -)+ g ?e/ 2 W . AM V b,(a;)^(r - r,) 

- w 2 - «7W ^ J w 2 , - w 2 - 47^ j 

Ne 2 /m , iNe/m , , .~ . . ,,. 

-» -5 5^— E(r,^) + — i__A c a;br,w) (44) 

in the approximation in which we assume the atoms are continuously distributed with a density iV and we 
define A c (u) = y/ p m h~fuj/ir, with p m — m/N . 
From Maxwell's equations and (|44p we obtain 

2 2 
V 2 E(r,w) + ^e(w)E(r,w) = -^K(r,w), (45) 



where the complex permittivity is defined as 



e(w) = 1 - -5 = f - -3 f-r— - e«(w) + ie/(w). 46 

We have also defined the "noise polarization" at frequency u>: 

K(r, W ) = AmNe l m AM b(r, W ). (47) 
— w z — 170-1 

This noise polarization obviously stems from the Langevin force Fx,j(i) in the quantum Langevin equation 
(|35p . Its principal properties for our purposes are the thermal equilibrium expectation values 



(K M (r,w)) = (^t(r, W ))=0, 
(/^(r,^(r',^')) = (^t(r,w)^(r , ,w'))=0, (48) 

and 

(i,t (r , w )^ (r >')> -4^ / H^(c,- W ')<5 3 (r-r') ? - 7 ^- T , (49) 



(£„(r, w)i^(r', w')> = 4n€j(w)^5(w - c/)5 3 (r - r') 



1 



a hui/k B T ^ 



(50) 



all of which follow from (j3"5)) and {bi^(uS)bj v {u)')) — 0. Equations (|4"9"| and (|ST))) are the well-known fluctuation- 
dissipation relations between the correlation functions of a noise polarization and the imaginary part of the 
dielectric function [TUl [H] . 

Next we define operators g\ (k, uj) by writing 

K(r,w) = fd 3 kY< 5A(k,^)e kA e lkr . (51) 

J A=l,2 

Since V ■ K(r, ui) — we can choose the vectors ekA such that k • ekA = 0, ekA ■ ekA' = 0, A = 1, 2; we also 
take the ekA to be real. Then 

g x (k,co) = [ rf 3 rK(r, W ).e k ,e- ik - r E (±-\ f d 3 rK x (v, u)e~^\ (52) 
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and Eqs. (|T7|) and Q28p imply the commutation relation 



We also define operators 
satisfying 



[g x (k,Lu),g{,(k',u;')] = ^ ei (u)S^5(u - Lo')S 3 (k - k'). (53) 

Z7T J 

c A (k,w) = [fte/M^vr 3 ]- 1 ^^) (54) 



[c A (k,w),4(k',o;')]=5AA^(w-a;')5 3 (k-k'). (55) 
Finally an expression for the quantized electric field follows from (14T1) . (|45|) . (|5T|) . and (|54|) : 

,2 /„2 



/• /■«> (.2 / 2 

E(r,t) = / d 3 fc V / ckjy/h^uj/2^- 2 / c A (k, o;)e k Ae- i ^ t - k - r ) + h.c. 

7 A -A) _ e(w)w 4 /<r 



(56) 



From V x E = — (l/c)<9B/9i we also obtain 

H(r, t)=i [d 3 kY] r duo^hej^)^- ^ _ . 2 c A (k, w) (k x e kA ) e -*(**-i<-') + h.c. (57) 

J A 7o « - Z\0J)0J Z /C 1 

These expressions have the same form as the corresponding ones obtained by Huttner and Barnett [TT] 
by Fano diagonalization of the entire system of coupled harmonic oscillators (EM field, dielectric oscillators, 
and bath oscillators). Their equations for the quantized electric and magnetic fields, however, involve 
annihilation and creation operators C* A (k, a;) and Ct(k,u)) for the exactly diagonalized Hamiltonian, instead 
of the reservoir annihilation and creation operators c A (k, oj) and ct (k, u>) appearing in our expressions (|56p 
and (|57[) . Their diagonalized Hamiltonian, including the zero-point energy, is 



H = I d 3 kY^ I dujhuj 



C5(k,w)C A (k,w) + l/2 . (58) 



The situation here parallels that for the simple model employed in Sections and [31 except that no rotating- 
wave approximation has been made, and that one deals with three coupled subsystems instead of two coupled 
subsystems: we again arrive at results by a straightforward "Langevin" approach that are equivalent to those 
obtained by diagonalizing the complete Hamiltonian. Equations (|56[) and (|5T[) are analogous to equation 
(flO|) obtained in the Langevin approach to the single oscillator coupled to a reservoir, whereas the Huttner- 
Barnett equations for the fields are analogous to equation (|25[) obtained by exact diagonalization. As in the 
model of Sections [2] and [3 results such as ([56]) and ([57]) obtained by the Langevin approach will reproduce 
those obtained by exact diagonalization for dielectric media in thermal equilibrium. To illustrate this we 
show in the next section that the total zero-point energy appearing in (|58p follows exactly from our approach; 
the calculation also sheds light on some of the physics involved, and in particular on the role of the Langevin 
forces in maintaining equilibrium. 



5 Energy Density 

We consider now the total energy density of the system of dielectric atoms, their reservoirs, and the electro- 
magnetic field, focusing for simplicity on the limit of zero temperature. We start from Poynting's theorem in 
the conventional notation, using the symmetrized Poynting operator S = (c/87r)[E xH-HxE] and taking 
expectation values over the initial state of the system consisting of the field, the dielectric atoms, and the 
reservoir: 

<S).„da = -^ /(E.^ + ^.E)^-J- /(H.^ + ^.H)^. (59) 



dt dt 1 8tt / v dt dt 
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The left-hand side gives the energy flux through a closed surface S and, given that we are assuming thermal 
equilibrium, must vanish. We identify the rate of change of the expectation value of the total energy density 
W as 

dW 1 . - <9D «9D Id . - ,. 

^^ (E -^ + ^' E) + 8^^ (H) ' (60) 
and the assumption of thermal equilibrium implies that this must also vanish. In the case of interest 
D = E + 47rP e + K, where P e is the part of the polarization giving rise to the dielectric permittivity e(uj) 
and K is the noise polarization defined by (jTT)) . Thus D = D e + K and 



dW _ dW x dW 2 

~dt ~ ~dT + ~bT' 

where 



and 



dW x 


= — (E 


<9D e 


<9D e 


•E) + 


1 d ~ 2 
8^dl {H 


dt 


dt 4 


dt 




dW 2 
dt 


— (E 


dK 


dK 

' ~dt 


E). 



(61) 

(62) 
(63) 



Before proceeding with the calculation of W we note the following identity expressing conservation of 
energy: 



dW d_ 

~dT ~ ^dt^ 



1 -2 1 

2 TOX i + 2 mw ° x J 



E 2 + H 2 



mk^-ij-FL^ir-rj)). (64) 



The first term is the rate of change of the energy density (kinetic plus potential) of the oscillators of the 
dielectric, while the second term is the rate of change of the energy density of the electromagnetic field. If 
there were no dissipation (7 = and therefore F^j = 0), the third term on the right would vanish, and W 
would be just the matter-plus-ficld energy density. The third term gives the rate of change of the energy 
density in the reservoirs; 2j^2j(^mic : j)5 3 (r — r^) is the rate of increase of energy density in the reservoir, 
equal to the rate at which the energy density of the dielectric oscillators decreases due to their coupling to 
their reservoirs, while ■ Xj • FLj5 3 (r — Yj) is the rate of work per unit volume done by the Langevin forces 
of the reservoirs on the dielectric oscillators. 
Using gTJ) and 

dt) r°° 

—^=-1 tto[e(w)E(r,w)e^ t -e*(uOE t (r,w)e^], (65) 
dt J 

and integrating over t, we obtain 

^M) = ^E f ^ - ^ (ftxfro,) • Et (r,J))e-^-^ + -1 (H 2 ( M )), (66) 

since the zero-temperature expectation value (E^(r, lo) ■ Fj\i (r, u>')) = and (E^(r, ui) ■ t\, (r, to')) = unless 
A = A' and u> = u>'. It is convenient to rewrite l|66p as a sum of two identical terms and to interchange u> and 
u) 1 in one of these terms; this allows us to write 

^M) = j~E /^^ " /£fl( i~^ ( "^ E A (r^).Et (r) a;0)e-^-^ 

OTT , Jq Jo IjJ — LO 

00 />oo 
du>' / did {lo' €l{ljj') + U!€l(u))) 



A 



8*^ 
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(E A (r, w) • e{(t, w'))e"* (a '" a '' ) * - (E A (r, c/) • E A (r, cj)^-"')' 1 .. ._, 



2(cj' - lu) 



8?r 



H'(r,t)>. 



(67) 



It follows from (|55[) and (|56p that the vacuum (zero temperature) expectation value 

(E A (r,c,)-Et(r,a;0) = (E A (r ) a;0-Ei(r,a;)) = A eiH ^ j ^___L__* (w _ ( 68 ) 

The first term in (|67p is now evaluated using 

- f .1 ^ r>(l,r\ f\ 

(69) 



hm — : — = — [we fl (w)]. 



We evaluate the second term by noting that the zeroth-order contributions in (lu — lu') in the numerator 
cancel each other, while the first-order terms result in a contribution linear in the elapsed time t: 



lim 



g— ui )t gi(cj— lj )t 



u 2(lu' — lu) 



it. 



(70) 



Therefore 

1 fi f°° ( d \ [ 



8tt 2 



\k 2 -£LU 2 /c 2 \ 2 87T 



^<H 2 (r,t)) 



(71) 



where we have used the relations e R = n 2 R — n 2 and ej = 2n R ni between the real and imaginary parts of the 
permittivity e(cu) and the refractive index n(u>). 

For the evaluation of W2 it is convenient to define K(k, lu) by writing 



K(r,t) = y du J d 3 fc^[K A (k,w)e-^ t e lk - r + K A (k,tj)e M *e- 4k - r ], 



and use (JSTJ), ([Ml), and dHU) to relate K A (k,w) and E A (k,< 



K A (k,w) = — [k 2 - eH W 7c 2 ]E A (k,^). 

LU Z 

After inserting (|72|) and ([73| in ([63]) and performing some algebra we get 



(72) 



(73) 



*- poo poo 2 I 

107T C Jg Jq CJ — W 



X / tffc 



-i(u—u')t 



i{u)—u')t 



k 2 — e(Lu)w 2 /c 2 k 2 — t* (lu)lu 2 / c 



and, proceeding as in the evaluation of Wi, 
W 2 (r,t) = - 



fOO poo L q a lo ' 

it) Re / duu' I du -J e/fcjWfct/Ww — to') 

c 2 ^ Jo Jo u- & 

/1 h f°° 



(74) 



(75) 
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We see that the time-dependent term in W%{t) exactly cancels the time-dependent term in W\(t). 
The total energy density is obtained by adding (fTTj) and ([75]) [15] : 



W 



8tt 2 c 3 



E 



dujuj s I Re 



aw 



a; (<w 



(76) 



Finally, using e(u;) = n 2 (cj) and the following relations 



Ree 3 / 2 



(n R - nj)n R + u;n R ( 2n R 
-- n 3 R - 3n R n 2 j, 



dm 



_ 2ni — 

du> duj 



ej T d 2 _ 2 (in/ 

— Im— (u ve) — 4n R nj + 2n R niuj— — , 
ui duj du) 

and summing over polarizations, we obtain the vacuum expectation value of the total energy density: 



W 



dn R \ 



2tt 2 c 3 



d 

du 1 



duJui 3 n R (ui) — [uin R (ijj)]. 



(77) 



(78) 



This has the exactly the form expected had we ignored absorption entirely and simply posited that each 
mode of frequency uj and wavenumber k = n R (uj)uj/c has a zero-point energy (l/2)?kj, so that the energy 
density summed over all modes is 



, i 

d k-Twj 



(79) 



The physical interpretation of this result is that the loss of energy due to absorption is balanced by the work 
done by the Langevin forces that maintain the canonical commutation relations of the dielectric oscillators. 



6 Summary 

Based on the simple model of a harmonic oscillator coupled to a reservoir, we showed how the Heisen- 
berg equations of motion leading to a Langevin equation for the oscillator can give results equivalent to 
those obtained from the exact (Fano) diagonalization of the complete oscillator-reservoir system. We then 
used the model of a dielectric medium as a collection of harmonic oscillators, each oscillator coupled to a 
reservoir responsible for dissipation and a Langevin force as well as to the electromagnetic field, to derive 
the fluctuation-dissipation relation between the noise polarization arising from the Langevin forces and the 
imaginary part of the permittivity of the dielectric medium. 

The simple oscillator-reservoir model we considered would suggest that the solutions for the electric and 
magnetic fields in a dielectric medium, with the noise polarization as a source, might have the same form 
as obtained when the complete system of dielectric oscillators, reservoirs, and the electromagnetic field is 
diagonalized. We showed that this is in fact the case. Then we considered the total energy density of the 
complete system and showed explicitly that it is given by Eq. ([75]) . which is exactly the form of the energy 
density obtained when absorption is ignored. In particular, we showed that a positive energy rate W\ > 
arising from the interaction of the electromagnetic field with the dielectric oscillators is exactly canceled 
by a corresponding negative energy rate coming from the interaction of the system with the reservoir, 
W 2 = -Wi < 0. 
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